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ABSTRACT 

The  acoustic  radiation  from  isotropic  turbulence  is  computed  numerically.  A  hybrid 
direct  numerical  simulation  approach  which  combines  direct  numerical  simulation  (DNS) 
of  the  turbulent  flow  with  the  Lighthill  acoustic  analogy  is  utilized.  It  is  demonstrated 
that  the  hybrid  DNS  method  is  a  feasible  approach  to  the  computation  of  sound  generated 
by  turbulent  flows.  The  acoustic  efficiency  in  the  simulation  of  isotropic  turbulence 
appears  to  be  substantially  less  than  that  in  subsonic  jet  experiments.  The  dominant 
frequency  of  the  computed  acoustic  pressure  is  found  to  be  somewhat  larger  than  the 
dominant  frequency  of  the  energy-containing  scales  of  motion.  The  acoustic  power  in  the 
simulations  is  proportional  to  eM®  where  c  is  the  turbulent  dissipation  rate  and  Mt  is  the 
turbulent  Mach  number.  This  is  in  agreement  with  the  analytical  result  of  Proudman 
(1952),  but  the  constant  of  proportionality  is  smaller  than  the  analytical  result.  Two 
different  methods  of  computing  the  acoustic  power  from  the  DNS  data  bases  yielded 
consistent  results. 

'This  research  was  supported  by  the  National  Aeronautics  and  Space  Administration  under  NASA 
Contract  No.  NASl-19480  while  the  authors  were  in  residence  at  the  Institute  for  Computer  Applications 
in  Science  and  Engineering  (ICASE),  NASA  Langley  Research  Center,  Hampton,  VA  24681. 
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Introduction 

In  1952,  Lighthill  posed  the  problem  of  estimating  the  sound  radiated  by  “a  fluctuating 
fluid  flow  occupying  a  limited  part  of  a  very  large  volume  of  fluid  of  which  the  remainder  is 
at  rest”  -  the  fluid  flow  of  course  is  either  ordered,  with  coherent  distributions  of  vorticity, 
or  disordered  and  turbulent.  It  is  only  now,  four  decades  later,  that  the  simplest  possible 
model  of  this  problem  i.e.,  sound  radiated  from  a  limited  body  of  isotropic  turbulence 
bounded  by  a  quiescent  fluid,  is  amenable  to  numerical  simulation. 

Isotropic  turbulence  in  a  bounded  domain  is  a  model  wherein  the  turbulence  is  unaf¬ 
fected  by  the  boundaries  enclosing  the  fluid,  and  furthermore  the  statistical  moments  are 
spatially  invariant  and  independent  of  orientation.  Isotropic  grid  turbulence  is  a  similar 
idealization,  in  that  the  turbulence  is  enclosed  by  wind  tunnel  walls  and  the  homogeneity 
of  the  turbulence  in  the  central  region  is  known  to  be  unaffected  by  the  wall  boundary 
layers.  The  problem  of  sound  emitted  by  isotropic  turbulence  has  been  investigated  by 
Proudman  (1952),  and  more  recently  by  Lilley  (1993)  in  an  analytical  study  and  by 
Witkowska  and  Juve  (1993)  in  a  large  eddy  simulation  (LES).  The  isotropic  turbulence 
is  represented  in  our  numerical  study  by  a  periodic  box  which  contains  the  fluid  motion 
and  is  surrounded  by  a  layer  in  which  the  velocities  decrease  to  zero.  The  fluid  motion 
inside  the  periodic  box  is  simulated,  and  it  is  assumed  that  the  sound  reuliated  by  the 
periodic  box  is  unaffected  by  the  surrounding  thin  layer. 

Two  possible  numerical  approaches  to  the  prediction  of  sound  produced  by  turbulence 
are  the  direct  numerical  simulation  (DNS)  and  the  hybrid  DNS.  In  a  DNS,  all  relevant 
scales  of  motion  are  numerically  resolved  in  a  solution  of  the  compressible  Navier-Stokes 
equations;  while  in  a  hybrid  DNS  only  the  turbulent  flow  is  resolved  numerically.  Sarkar 
and  Hussaini  (1993)  have  estimated  the  turbulence  Reynolds  numbers  Ret  that  could 
be  obtiuned  for  a  three-dimensional  calculation  of  sound  generated  by  turbulence  with 
the  two  methods  for  a  given  spatial  grid.  It  was  found  that  for  Aft  =  0.01  and  0.1,  the 
Ret  that  can  be  simulated  by  DNS  of  sound  on  a  uniform  grid  is  respectively  0.004  and 
0.08  times  smaller  than  that  achieved  in  the  hybrid  DNS  method.  Evidently,  the  much 
larger  computational  volume  required  in  the  DNS  of  sound  approach  leads  to  a  much 
lower  turbulence  Reynolds  number  relative  to  the  hybrid  approach.  By  stretching  the 
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computational  grid  or  by  using  overlapping  grids,  the  Ret  obtainable  with  the  DNS  of 
sound  approach  can  be  increased.  However,  the  effect  of  grid  stretching  on  the  accuracy 
of  the  far-field  acoustic  pressure  which  is  several  orders  of  magnitude  smaller  than  the 
aerodynamic  pressure  has  to  be  carefully  examined.  Crighton  (1993)  points  out  that 
computations  of  low  Mach  number  aeroacoustics  by  the  DNS  of  sound  approach  have 
to  contend  with  a  number  of  difficulties  which  include  the  large  disparity  between  the 
acoustic  length  scale  and  the  turbulence  integral  length  scale,  the  much  smaller  energy  in 
the  acoustic  far  field  relative  to  the  turbulent  energy,  and  the  multipole  structure  of  the 
basic  acoustic  source.  Due  to  these  reasons,  Lighthill  (1993)  recommends  the  use  of  an 
acoustic  analogy  for  the  computation  of  sound  radiated  from  low  Mach  number  flows.  The 
present  work  employs  the  the  hybrid  DNS  approach  with  the  Lighthill  acoustic  analogy 
for  the  purpose  of  computing  the  sound  generated  by  three-dimensional,  fully  turbulent 
flow.  Lighthill’s  analogy  has  been  used  previously  to  compute  the  sound  numerically  from 
the  large-scale  instability  waves  in  a  shear  layer  by  Gatski  (1979),  but  we  are  not  aware 
of  any  previous  attempt  to  use  the  analogy  within  the  framework  of  three-dimensional 
DNS  of  turbulence. 

The  2K;oustic  analogy  which  was  first  proposed  by  Lighthill  (1952)  reduces  the  aeroa¬ 
coustics  problem  to  an  inhomogeneous  wave  equation.  In  his  analogy,  Lighthill  chose  the 
homogeneous  part  to  be  a  linear  wave  equation  for  the  density  and  obt^ned  the  acoustic 
far-held  as  an  integral  of  the  equivalent  acoustic  sources  in  a  uniform  medium  at  rest 
which  replace  the  fluid  motion  inside  a  bounded  volume.  Furthermore,  by  considering 
the  statistics  of  the  acoustic  sources,  Lighthill  deduced  the  eighth-power  law  at  low  Mach 
numbers  for  the  net  acoustic  power  radiated  from  a  turbulent  region 

Lighthill’s  analogy  was  followed  by  other  analogies  (Phillips  (1960),  Powell  (1964), 
Ribner  (1964),  Lilley  (1974),  and  Howe  (1975)).  In  the  analogies  of  Phillips,  Lilley  and 
Howe,  the  homogeneous  part  of  the  equation  for  the  acoustic  variable  is  a  nonlinear 
convected  wave  equation  rather  than  the  linear  wave  equation  in  Lighthill’s  analogy, 
and  as  discussed  by  Goldstein  (1974),  Crighton  (1975)  and  Ffowcs  Williams  (1977)  has 
the  advantage  of  explicit  consideration  of  the  convection  and  refraction  of  sound  by 
turbulence.  Powell  and  Howe  rewrite  the  term  u-  Vu  as  u;  x  u-|-  Vu*/2  in  their  analogies 
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which  leads  to  an  explicit  vorticity-dependent  term  in  the  inhomogeneous  part  of  the  wave 
equation.  The  explicit  appearance  of  the  vorticity  in  the  forcing  term  is  convenient  in 
unsteady  flows  such  as  the  Von  Karman  vortex  shedding  in  the  laminar  wake  of  a  cylinder, 
but  perhaps  not  so  important  in  high  Reynolds  number  turbulent  flows  where  the  flow  is 
more  incoherent  and  cannot  be  idealized  as  a  simple  collection  of  idealized  vortices.  It  is 
also  more  convenient  in  the  context  of  vortex  methods  (Leonard  (1985)).  Applications  of 
the  Lighthill  theory  are  consistent  with  experimental  data  on  noise  radiated  from  subsonic 
jets  running  at  ambient  temperatures.  However  in  supersonic  jets  (Seiner(1992)),  it 
may  be  necessary  to  account  for  additional  physical  features  explicitly  such  as  Mach 
wave  radiation,  jet  screech,  and  other  shock-associated  noise.  In  our  work  on  noise 
generated  from  isotropic  turbulence  at  low  Mach  numbers,  Lighthill’s  acoustic  analogy 
is  the  preferred  representation  and  is  the  one  followed  herein. 


Computation  of  the  Lighthill  acoustic  analogy 


The  acoustic  analogy  of  Lighthill  can  be  written  in  the  space  -  derivative  form 


r 


dy, 


(1) 


but  for  a  far-field  location  of  the  observer  point  with  |x|  >>  \y\,  Lighthill  showed  that 
the  integrand  in  (1)  can  be  approximated  by  the  second  time  derivative  taken  at  the 
retarded  time  yielding 

-  r/co)dy,  (2) 

where  pa  is  the  fluctuation  of  pressure  relative  to  the  ambient,  co  is  the  ambient  speed 
of  sound,  and  r  =  |x  —  y|.  In  Eq.  (1),  density  fluctuations  are  neglected  in  the  source 
term  and  the  approximation,  ~  d^{puiUj) I dyidyj.,  is  utilized  since  only  low  Mach 
numbers  are  considered  in  the  present  work.  Here,  p  is  the  mean  density  which  is  of 
course  a  constant  equal  to  the  ambient  density  of  the  surrounding  medium  which  is  at 
rest.  In  Eq.  (2),  Tij{y,t~rlco)  denotes  d^{pUiUj)ldt‘^  evaluated  at  retarded  i\met  —  r/cQ. 

The  major  advantage  of  the  time-derivative  form  of  the  acoustic  analogy  over  the 
space-derivative  form  is  that  the  quadrupole  nature  of  the  sound  source  is  directly  evident 
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in  Eq.  (2)  but  not  in  Ek|.  (1).  However,  the  computer  storage  requirements  of  the  space- 
derivative  form  is  smaller  by  a  factor  of  6  relative  to  those  of  the  time-derivative  form 
which  could  present  a  significant  advantage  in  large-scale  computations  where  computer 
memory  is  a  constraint.  In  Sarkar  and  Hussaini  (1993),  the  relative  axlvantages  of  Ek}s.  (1) 
and  (2)  were  investigated  for  a  test  problem  with  an  imposed  quadrupole.  The  retarded 
time  effect  was  accounted  for  by  the  time  accumulation  method  wherein  the  observer 
time  t  was  approximated  by  <  ~  [(r  -f  cor)/At]At  where  [.]  denotes  the  integer  value 
function,  and  At  denotes  the  time  step  used  for  time  advancement  of  the  flow.  It  was 
found  that  the  number  of  time  points  required  for  the  spatial  derivative  form  is  too  large; 
turbulence  of  the  quadrupole  source  type  requires  0{lfMf)  time  points  per  oscillation 
of  the  source.  Due  to  the  stringency  of  this  requirement,  we  choose  the  time-derivative 
form,  Eq.  (2),  over  the  space-derivative  form,  Eq.  (1). 

Furthermore,  Eq.  (1)  is  not  quite  appropriate  here  because  of  the  periodic  boundary 
conditions  used  for  the  turbulence  simulation.  Because  of  the  non-zero  velocity  at  the 
fluid  boundary  of  the  periodic  box,  the  dipole  contribution  to  the  sound  from  the  surface 
of  the  flow  domain  is  non-zero  and  would  be  included  in  a  calculation  using  Ekj.  (1). 
The  surface  contribution  is  an  artifact  of  the  boundary  conditions  and  is  not  of  interest 
to  us.  Eq.  (2),  on  the  other  hand,  gives  directly  the  quadrupole  sound  characteristic 
of  a  finite  volume  of  turbulent  flow  embedded  in  a  medium  at  rest.  We  ensure  that  the 
periodic  boundary  conditions  do  not  unrealistically  affect  both  the  flow  field  and  acoustic 
source  by  keeping  the  length  of  the  computational  domain  much  larger  than  the  spatial 
correlation  lengths  of  the  velocity  field  Uj  and  the  acoustic  source  Tij. 

In  the  case  of  unforced  isotropic  turbulence,  the  acoustic  source  term,  /  T',j(y,  t  —  r/co)  dy, 
decays  in  time.  Consequently,  the  acoustic  pressure  is  statistically  unsteady  and  time  av¬ 
eraging  is  inappropriate  for  statistical  analysis.  Short-time  averaging  was  used  in  Sarkar 
and  Hussaini  (1993)  to  compute  the  statistics  of  the  acoustic  pressure.  But  such  a  method 
has  since  proven  to  be  unsatisfactory  due  to  the  time  scale  of  the  turbulence  decay  being 
of  the  same  order  as  the  time  scale  of  the  acoustic  fluctuations.  The  results  presented 
here  are  substantially  different  from  the  preliminary  work  of  Sarkar  and  Hussaini  (1993). 

In  the  present  work,  ensemble  averages  are  used  to  compute  the  acoustic  intensity,  power 
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and  frequency  spectrum.  20  cases  were  run  with  the  same  parameters  and  the  same 
initial  statistics,  but  with  different  instantaneous  initial  conditions.  For  each  case,  the 
observation  points  were  distributed  on  two  concentric  spheres  in  the  far-field,  10  on  each 
sphere.  Due  to  the  isotropy  of  the  turbulence,  the  20  different  simulations  with  10  dif¬ 
ferent  observer  points  at  the  same  distance  from  the  center  of  the  flow  domain  lead  to 
200  different  samples  of  the  acoustic  pressure  which  is  a  sufficiently  large  number  to  ob¬ 
tain  acceptable  statistics.  The  scaled  au^oustic  power  is  obtained  by  averaging  over  400 
samples. 

Flow  simulation  method 

The  turbulent  flow  inside  a  cubical  domain  is  obtained  by  solving  the  compressible 
Navier-Stokes  equations  by  a  numerical  algorithm  which  was  originally  developed  to  in¬ 
vestigate  compressibility  effects  in  isotropic  turbulence  (Erlebacher,  Hussaini,  Kreiss,  and 
Sarkar  (1990),  Sarkar,  Erlebacher,  Hussaini,  and  Kreiss  (1991a)),  and  homogeneous  shear 
turbulence  (Sarkar,  Erlebacher,  and  Hussaini(1991b),  Sarkar  (1992)).  The  homogeneous 
turbulence  problem  permits  periodic  boundary  conditions  in  all  three  coordinate  direc¬ 
tions.  It  is  of  course  necessary  that  the  length  of  the  computational  domain  be  much 
larger  than  the  integral  length  scale  of  the  turbulence  for  obtaining  realistic  flow  fields. 
Spectral  accuracy  is  obtained  by  using  a  Fourier  collocation  method  for  the  spatial  dis¬ 
cretization  of  the  governing  equations.  A  third  order,  low  storage  Runge-Kutta  scheme 
is  used  for  advancing  the  solution  in  time. 

Initial  conditions  have  to  be  prescribed  for  u/,  p,  p  and  T.  The  initial  velocity  field 
is  split  into  two  independent  components,  that  is,  u/  =  uj'  uf\  where  each  component 
has  a  zero  average.  The  solenoidal  velocity  field  uf  which  satisfies  =  0  is  chosen 

to  be  a  random  Gaussian  field  with  the  power  spectrum 

E{k)  =  k*  exp{—2k^ /kl^)  (3) 

where  km  denotes  the  wave  number  corresponding  to  the  peak  of  the  power  spectrum. 
The  compressible  velocity  uf  which  satisfies  V  x  =  0  is  also  chosen  to  be  a  random 
Gaussian  field  satisfying  the  same  power  spectrum,  Eq.  (3).  The  power  spectra  of  the 
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two  velocity  components  are  scaled  so  as  to  obtain  a  prescribed  Umu  =  and 

a  prescribed  x  =  «^/«nn*  which  is  the  compressible  fraction  of  kinetic  energy.  The 
pressure  p^'  associated  with  the  incompressible  velocity  is  evaluated  from  the  Poisson 
equation 

(4) 


=  —pu 


V  V 
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It  remains  to  specify  the  initial  values  of  the  thermodynamic  variables.  The  mean  density 
p  is  chosen  equal  to  unity,  and  p  is  chosen  so  as  to  obtain  a  prescribed  Mach  number 
«nns/Y7P/^  characterizing  the  turbulence.  The  fluctuating  density  p'  and  compressible 
pressure  are  chosen  as  random  fields  with  the  power  spectrum  determined  by  Eq.  (3) 
and  prescribed  pms  and  p^.  The  pressure  then  becomes  p  =  p  +  p^'  +  p^',  the  density 
is  p  =  p  +  p',  and  the  instantaneous  temperature  T  is  obtained  from  the  equation  of 
state  p  =  pBT.  We  note  that  the  splitting  of  velocity  and  pressure  into  incompressible 
and  compressible  components  is  confined  to  the  specification  of  initial  conditions.  No 
such  split  is  carried  out  while  computing  the  acoustic  sources  during  the  time  evolu¬ 
tion.  Although  the  numerical  algorithm  is  capable  of  simulating  compressible  turbulence 
with  non-zero  density  fluctuations  and  non-zero  dilatation,  we  consider  low  Mach  num¬ 
ber,  qucisi-incompressible  turbulence  in  the  present  paper.  The  €u:oustic  source  Tij  is 
approximated  by  puiUj  and  density  fluctuations  neglected  in  the  acoustic  source. 


Characteristics  of  the  simulated  turbulence 

The  aicoustic  radiation  from  isotropic  turbulence  has  been  computed  with  the  hybrid 
DNS  method  for  the  case  with  the  following  initial  parameters 

Rex  =  50,  ,  Mt  =  0.05 

where  Rex  is  the  Taylor  microscale  Reynolds  number,  and  Mt  the  turbulent  Mach  num- 
ber.(  Note  that  Rex  =  <l^lv  where  q  =  yju\u\,  X  =  is  the  fluctuating 

vorticity,  and  v  is  the  kinematic  viscosity;  while  Mt  =  q/c  where  c  is  the  mean  speed 
of  sound.)  The  initial  energy  spectrum  is  given  by  Eq.  (3)  with  km  =  6.  The  initial 
data  for  the  case  discussed  here  is  chosen  to  be  incompressible;  the  velocity  is  solenoidal 
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(x  =  0),  the  density  is  constant,  and  the  pressure  is  initialized  with  the  usual  Poisson 
equation  applicable  to  incompressible  flows.  Twenty  simulations  on  a  64^  grid  with  a 
time  step  of  At  =  0.00275(A'/c)o  were  caried  out.  A  single,  higher- resolution  simulation 
was  performed  on  a  128^  spatial  grid  with  At  =  0.001375(^/6)0,  other  parameters  re¬ 
maining  equal,  to  check  the  accuracy  of  the  coarse-grid  solution.  The  simulations  were 
conducted  for  one  eddy  turnover  time  (eot/A'o  =  1  )•  The  64^  simulations  were  used 
to  obtain  ensemble- averaged  statistics  from  the  instantaneous  acoustic  pressure,  while 
the  128^  simulation  was  used  to  obtain  the  acoustic  power  from  the  appropriate  flow 
statistics. 

The  turbulence  statistics  at  a  given  time  in  a  DNS  case  are  obtained  by  averaging  over 
the  computational  volume.  For  example  the  turbulent  kinetic  energy  K  =<  u,u,  >  /2 
where  <  .  >  denotes  a  volume  average.  In  the  present  simulations,  the  turbulence  is 
unforced.  Consequently,  the  turbulent  kinetic  energy  decays  in  time.  Figs.  1  and  2  show 
the  evolution  of  the  turbulent  kinetic  energy  K  and  turbulent  dissipation  rate  e,  respec¬ 
tively.  K  decreases  by  about  a  factor  of  5  by  the  end  of  the  simulation.  The  evolution  of 
K  for  tot/Kf,  >  0.2  is  well  represented  by  a  power  law  decay  (t  —  <o)~”  with  the  exponent 
n  =  1 .4,  a  value  which  is  slightly  larger  than  experimentally  measured  exponents  which 
lie  in  the  range  1.15  <  n  <  1.35.  The  turbulent  dissipation  rate  e  increases  during  the 
initial  transient  in  Fig.  2  due  to  the  generation  of  small-scale  fluctuations,  and  finally 
decays  in  the  absence  of  external  forcing.  The  microscale  Reynolds  number  Rex  whose 
evolution  is  shown  in  Fig.  3  decays  by  a  factor  of  3  at  the  end  of  the  simulation.  The 
simulation  is  terminated  at  a  nondimensional  time  eotfKo  =  1,  because  Rex  at  later 
times  could  become  too  small  to  represent  realistically  the  spatio-temporal  complexity 
of  turbulence.  The  skewness  of  the  velocity  derivative  Sk  =  (9u/9i)^/[(5u/dx)^]^/^  is 
a  measure  of  the  non-linear  vortex  stretching.  According  to  Lesieur(1990),  experiments 
on  grid  turbulence  give  Sk  —0.4,  while  simulations  give  Sk  ~  —0.5.  Fig.  4  shows  that 
skewness  factor  asymptotes  to  a  value  of  Sk  =  —0.46  in  the  DNS  which  is  consistent 
with  previous  experiments  and  simulations. 
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Characteristics  of  the  computed  acoustic  pressure 

The  fluctuating  acoustic  pressure  is  computed  from  the  DNS  data  by  the  temporal 
form  of  the  Lighthill  analogy, 

where  [.]  denotes  evaluation  at  retarded  time,  while  the  overbar  denotes  ensemble  average. 
Thus,  pa{x.,  t)  is  computed  from  Eq.  (2);  it’s  ensemble  average  calculated  over  the  20  DNS 
Ccises;  and  finally,  the  ensemble  average  is  subtracted  out  to  calculate  p^(x,  t).  Fig.  5 
shows  the  evolution  of  the  acoustic  pressure  fluctuation  normalized  by  the  ambient 
pressure  Pq.  The  smallness  of  the  acoustic  pressure  fluctuation  which  is  about  seven 
orders  of  magnitude  smaller  than  the  thermodynamic  pressure  would  present  a  formidable 
challenge  for  the  direct  computation  of  sound  by  simulation  of  the  compressible  Navier- 
Stokes  equations  in  the  far-field.  Since  the  turbulent  source  decays  with  time,  the  acoustic 
pressure  also  decays  with  time.  We  have  tried  to  factor  out  the  decay  in  turbulence 
statistics  by  dividing  the  random  pressure  signal  p^(t)  by  an  appropriate  deterministic 
function  g{t).  The  choice  of  g{t)  was  guided  by  the  analysis  of  Proudman  (1952)  who 
considered  the  generation  of  noise  by  isotropic  turbulence  and  used  statistical  models  of 
various  two-point  moments  within  the  framework  of  the  Lighthill  analogy  to  obtain  the 
following  expression  for  the  acoustic  power  Pa  per  unit  mass 

Pa  =  a(«V0^  (6) 

where  a  is  a  constant  related  to  the  shape  of  the  longitudinal  velocity  correlation  /(r),  u 
is  the  root  mean  square  (rms)  of  o^t  of  the  velocity  components,  and  /  is  the  longitudinal 
integral  length  scale  of  the  velocity.  For  the  Heisenberg  form  of  the  energy  spectrum, 
Proudman  obtains 

Pa  =  38(«V0^  (7) 

% 

and,  with  the  assumption  /(r)  =  Proudman  obtains 

Pa  =  13(«V/)^  (8) 
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Eq.  (8)  implies  that  g{t)  cx  yjtMf  can  be  used  to  normalize  the  acoustic  pressure  pj^ 
in  order  to  factor  out  the  temporal  decay.  Fig.  6  shows  the  evolution  of  rescaled  <icoustic 
pressure  p\  defined  by 

*  ^{€LI4)m} 

Here,  L  is  the  length  of  the  computational  region.  It  appears  that  the  normalization 
does  factor  out  the  temporal  decay  of  the  acoustic  pressure.  Fig.  7  shows  the  acoustic 
pressure  fluctuations  at  different  observer  points.  The  instantaneous  pressure  signal  at  a 
given  observation  point  is  non-periodic,  chaotic,  and  has  a  range  of  frequencies. 

The  rms  acoustic  pressure  Prms(0  observed  at  a  distance  x  is  obtained  by  ensemble¬ 
averaging  over  200  samples  of  fluctuating  acoustic  pressure  These  samples  are  ob¬ 
tained  from  the  20  cases  in  the  DNS  with  each  case  having  10  observation  points  per 
sphere  distributed  on  two  concentric  spheres  around  the  center  of  the  flow  domain.  The 
acoustic  power  emitted  from  the  volume  V  of  fluid  is 

Pv{t)  = 

PqCo 

and  the  a^^oustic  power  per  unit  mass  of  turbulent  fluid  is 

/ox 


PaH)  = 


PoV  • 


Fig.  8  shows  the  acoustic  power  in  decibels  from  the  DNS  and  the  theoretical  result. 
The  dotted  curve  in  Fig.  8  corresponds  to  a  value  of  a  =  13  in  Proudman’s  work  giving 

Pa  =  8.7£(i)^ 

=  0.5e(f)Mt(t)®  (10) 

and  is  obtained  from  Eq.  (8)  by  using  the  definition  of  Mt  =  y/2K/co  and  our  DNS 
result  e  ~  l.5u^fl.  The  computed  power  is  smaller  than  that  obtained  from  Eq.  (10)  by 
6  —  8  dB  over  the  time  of  the  simulation.  The  ratio  FA/eMf  varies  between  0.15  and 
0.07  in  Fig.  9,  i.e.,  by  a  factor  of  2  during  the  decay  of  turbulent  kinetic  energy,  while 
the  numerator  and  denominator  vary  by  about  a  factor  of  15  during  that  time.  Thus, 
the  DNS  is  consistent  with  Proudman’s  result  that  Pa  is  proportional  to  eMf .  However, 
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if  the  proportionality  constant  from  the  DNS  is  taken  to  be  ~  0.10,  it  is  smaller 

relative  to  Proudman’s  result,  Eq.  (10),  by  a  factor  of  5. 

Witkowska  and  Juve  (1993)  have  computed  the  sound  generated  from  unforced  isotropic 
turbulence  by  performing  large  eddy  simulations  (LES)  on  16^  spatial  grids  for  determin¬ 
ing  the  turbulent  flow  and  then  using  the  Lighthill  acoustic  analogy  for  computing  the 
far-field  sound.  Five  turbulent  Mach  numbers  in  the  range  of  0.012  to  0.015  were  con¬ 
sidered;  the  initial  Rex  was  about  400.  In  agreement  with  our  DNS  results,  the  LES 
study  finds  that,  although  the  acoustic  power  is  proportional  to  eA/®,  it  is  12  dB  less 
than  Proudman’s  result  for  the  Heisenberg  spectrum.  According  to  Witkowska  and  Juve 
(1993),  a  reconsideration  of  Proudman’s  analysis  with  the  von  Karman  spectrum  instead 
of  the  Heisenberg  spectrum  decreases  the  theoretical  estimate  so  that  the  LES  result  is 
6dB  lower  than  the  revised  theoretical  result. 

The  turbulence  in  the  DNS  has  moderate  to  low  Reynolds  number  since  Rex  decays 
from  a  value  of  50  to  15  during  the  simulations.  The  LES  of  Witkowska  and  Juve  (1993) 
has  a  turbulence  Reynolds  number  which  is  a  factor  of  10  higher  than  the  DNS,  albeit  with 
the  approximations  inherent  in  a  subgrid  scale  model.  Although,  there  is  a  difference  in 
Reynolds  number,  the  DNS  and  LES  studies  are  consistent  in  giving  an  acoustic  power 
which  is  a  factor  of  5  smaller  than  Proudman’s  result  Eq.  (10).  Proudman’s  analysis 
predates  the  present  computational  studies  by  about  four  decades  and,  of  necessity,  had 
to  employ  vcirious  assumptions  such  as  Gaussian  statistics  for  the  velocity  and  its  time 
derivatives,  neglect  of  retarded  time  during  evaluation  of  the  integral  expression  for  the 
acoustic  power,  neglect  of  the  decay  of  unforced  isotropic  turbulence,  and  a  specific 
shape  for  the  longitudinal  space  correlation  /(r),  Lilley  (1993),  in  a  reconsideration  of 
Proudman’s  analysis,  finds  that  the  coefficient  a  in  Eq.  (6)  is  likely  to  be  between  3  and 
10.  In  contrast,  our  DNS  gives  a  ~  2.6.  The  analysis  of  Lilley  (1993)  also  elegantly 
shows  that  the  acoustic  power  is  sensitive  to  the  shape  of  the  space-time  covariance  of 
Txx,  the  flatness  factor  of  the  velocity,  and  is  proportional  to  the  fourth  power  of  the 
characteristic  Strouhal  number  s  =  ul/u  where  /  and  u  are  the  integral  scales  associated 
with  the  spatial  and  temporal  autocorrelations  of  the  velocity,  respectively,  and  u  is  the 
rms  of  a  velocity  component.  Using  the  velocity  flatness  factor  of  3.0  in  our  DNS  results. 
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Lilley  (1993)  finds  a  =  3.4  with  5  =  1  and  a  —  8.3  with  5  =  1.25  in  an  analytical 
calculation  which  assumes  the  same  f{r/l)  that  corresponds  to  Proudman’s  a  =  13. 
Given  the  sensitvity  of  a  to  the  assumed  model  correlations,  it  is  not  surprising  that 
the  analytically  determined  values  for  or  are  somewhat  different  from  the  corresponding 
results  in  the  simulations.  Overall,  the  DNS  value  of  a  is  consistent  with  the  lower  end 
of  the  range  of  values  of  a  found  by  Lilley  (1993). 

The  scaled  acoustic  power  is  significantly  smaller  than  that  typical  of  sub¬ 

sonic  jets.  In  subsonic  jets  the  acoustic  eflSciency  (sound  power/jet  power)  according  to 
Lighthill  (1954)  is  approximately  10“"^ Af®  ,  where  M  is  the  jet  Mach  number.  The  tur¬ 
bulence  intensity  UrmsIUjet  can  be  roughly  estimated  to  be  0.17,  and  c  is  roughly  0.1  of 
the  jet  power,  which  implies  that  Pyi/cA/f  ~  7.  Thus  the  sound  radiated  from  isotropic, 
homogeneous  turbulence  is  significantly  smaller  than  that  typical  of  jet  acoustics.  We 
argue  below  that  one  of  the  reasons  for  this  dissimilarity  is  the  homogeneity  constraint. 

The  acoustic  pressure  fluctuation  is  given  by  Eq.(5).  Consider  the  integral  in  Eq.  (5) 
as  Mt  0. 

J  Tij  -  fij  dy 

=  V  <  fij  >  -V<  fi,  >.  (11) 

The  second  line  in  Eq.  (11)  follows  because  retarded  time  can  be  neglected  in  the  limit  of 
M(  — >  0,  and  the  third  line  follows  from  the  definition  of  the  volume  average  denoted  by  < 

.  >.  In  the  case  of  homogeneous  turbulence,  <  Tij  >  is  equivalent  to  an  ensemble  average, 
and  since  further  time  averaging  does  not  change  a  quantity  which  is  already  an  average 

<  Tij  >  =<  Tij  >.  Consequently  the  integral  in  Eq.  (5)  is  zero  when  retarded  time  effects 
are  neglected  in  the  limit  of  Mt  — >  0.  However,  in  the  case  of  inhomogeneous  turbulence 

<  Tij  >  is  not  a  temporally  smooth  quantity,  and  this  integral  is  not  necessarily  zero 
even  if  retarded  time  effects  were  neglected. 

For  small  but  finite  Mach  number,,  when  retarded  time  effects  cannot  be  neglected, 
it  is  clear  that  the  volume  integration  smoothes  out  the  temporal  oscillation  in  [T.j]  to  a 
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greater  extent  in  homogeneous  turbulence  than  in  inhomogeneous  turbulence,  which  thus 
causes  smaller  fluctuations  in  the  radiated  sound  from  homogeneous  turbulence  relative 
to  the  inhomogeneous  case. 

The  frequency  content  of  the  acoustic  pressure  was  obtained  by  analyzing  the  nor¬ 
malized  acoustic  pressure  p\{t)  in  Fourier  space.  The  power  spectrum  Ep{ij)  obtained  by 
averaging  over  400  records  of  p\{t)  is  shown  in  Fig.  10.  The  frequency  is  normalized  by 
/o/uo>  where  Iq  is  the  initial  longitudinal  integral  length  scale  and  uo  is  the  initial  rms  of 
a  velocity  component.  The  power  spectrum  in  Fig.  10  has  a  peak  at  a  Strouhal  number 
St  =  ujlofuo  ~  3.5.  Thus,  the  dominant  frequency  of  the  sound  radiated  to  the  far-field 
in  the  isotropic  turbulence  DNS  is  somewhat  larger  than  the  dominant  frequency  uq/Iq 
of  the  energy-containing  range  of  turbulence. 

Computation  of  the  acoustic  power  from  the  turbulence  statistics 

The  acoustic  power  obtained  in  our  simulations  is  somewhat  smaller  than  predicted 
by  Proudman’s  analysis.  This  discrepancy  motivated  an  alternative  method  to  obtain 
the  cicoustic  power  wherein  the  the  fourth-order,  two-point  correlation  T^i(0)T^j.(r)  from 
the  DNS  is  used. 

Consider  an  observer  point  at  location  (x,0,0).  Then,  due  to  Eq.  (5)  for  the  acoustic 
pressure  fluctuation,  the  acoustic  power  Pa  per  unit  mass  at  the  observer  point  given  by 
Eq.  (9)  becomes 

PA{t)  =  j  j  dx^dxB  (12) 

where  V  is  the  source  volume,  and  x^  and  tA  denote  the  position  of  source  point  A  and 
time  of  sound  emission  at  point  A  so  as  to  reach  observer  at  time  t;  similarly,  Xb  and  Ib 
refer  to  source  point  B.  Proudman  assumes  that  7’^^(x^,i)7’^jp(xB,t)  is  constant  over  the 
maximum  time  difference  tA  —  te  so  that  retarded  time  effects  can  be  ignored  during  the 
volume  integration  with  respect  to  source  point  B.  This  assumption  leads  to 

J  f^^i-KA,tA)TLi^BjA)  dXA  .  (13) 
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Proudman,  then  assumes  that  can  be  considered  constant  during  the  integration  over 
source  point  A,  to  obtain  finally 

PaW  =  ^  I  fU^A)fi.M  dXA  (14) 

The  two-point  spatial  autocorrelation  in  Elq.  (14)  is  a  function  of  r  =  |x^— XbI  only,  due  to 
the  turbulence  isotropy.  Let  us  define  aj  =  («*—  <  <  u*  >,w^—  <  >),  and 

denote  the  spatial  correlation  by  Qy(r).  The  subscript  t  =  1,2,3 

denotes  the  x,y,z  components,  respectively.  Then,  Eq.  (14)  becomes 

PaW  =  -^j  <?n(r)*  (15) 

The  following  expression  for  Qij  can  be  written  due  to  the  isotropy  constraint  (e.g.  see 
Hinze(1975)) 

g«(r)  =  +  G(r)<«l 

where  the  fluctuation  {d^u^jdt^y  =  d^u^fdt^—  <  d^u^jdt^  >,  while  F{r)  and  G{r)  are 
the  longitudinal  and  lateral  correlation  functions  associated  with  Qij  respectively.  Thus, 

Qii(r)  =  {dhJ/dt^  (F(r)  cos*(0)  +  G{r)  sin*(0)] 


where  d  is  the  angle  between  r  and  the  x  axis.  The  volume  integration  in  Eq.  (15)  with 
respect  to  the  separation  vector  r  is  then  performed  using  a  spherical  coordinate  system 
to  yield  _ 


PA{t)  = 


(d^uydt^) 

4 


2V2  ,«_2 


/„  j('"('')  +  2G(r))*  , 


(16) 


Proudman  (1952)  uses  the  statistical  theory  of  isotropic  turbulence  to  model  the  corre¬ 
lations  in  Ekj.  (16);  we,  on  the  other  hand,  directly  compute  F{r)  and  G(r) 

from  the  DNS  data  base. 

The  longitudinal  correlation  functions  F(r)  associated  with  Qij,  and  /(r)  associated 
with  the  velocity  field  are  shown  in  Fig.  11  at  the  time  Cot/Fo  =  0.7.  Evidently,  the 
presence  of  the  second  time  derivative  leads  to  a  much  smaller  correlation  length  of 
the  acoustic  source  relative  to  the  velocity.  To  confirm  our  results,  a  new  simulation 
was  performed  with  twice  the  resolution  in  both  space  and  time.  Fig.  12  shows  that  the 
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higher  resolution  128^  simulation  gives  a  result  for  F(r)  which  is  similar  to  the  coarse-grid 
simulation.  The  correlation  length  scale  for  the  acoustic  source  is  0.10,  while  the  integral, 
Taylor  and  Kolmogorov  length  scales  are  0.58, 0.12,  and  0.03,  respectively.  It  thus  appears 
that  the  correlation  length  of  in  our  simulation  of  isotropic  turbulence  is  smaller 

than  the  turbulence  integral  length  scale  and  larger  than  the  dissipative  length  scale. 

Eq.  (16)  W21S  used  to  provide  an  alternative  estimate  of  the  acoustic  power  from  the 
DNS.  Fig.  13  shows  that  the  result  from  Ek|.  (16)  is  consistent  with  the  previous  result 
obtained  from  statistical  post-processing  of  the  instantaneous  acoustic  pressure.  Since 
retarded  times  are  neglected  in  simplifying  the  exact  expression  for  the  acoustic  power 
Eq.(12)  to  obtain  Eq.  (16),  it  is  not  surprising  that  there  is  some  difference  between  the 
two  methods. 

Conclusions 

We  have  applied  the  hybrid  DNS  approach  to  the  problem  of  sound  generated  by 
isotropic  turbulence.  The  combination  of  a  spectral  DNS  of  isotropic  turbulence  and 
the  Lighthill  acoustic  analogy  is  demonstrated  to  be  a  feasible  approach  for  obtaining 
the  instantaneous  acoustic  pressure.  Although  the  time  derivative  form  of  the  Lighthill 
acoustic  analogy  requires  more  computer  memory  for  storage  than  the  spatial  derivative 
form,  it  is  preferable  due  to  a  less  stringent  time  step  constraint. 

The  numerically  computed  acoustic  efficiency  of  isotropic  turbulence  is  much  smaller 
than  that  in  subsonic  jet  experiments.  The  frequency  spectrum  of  the  acoustic  pressure 
indicates  that  the  dominant  acoustic  frequency  is  somewhat  larger  than  the  frequency 
of  the  energy-containing  eddies  as  found  in  Proudman  (1952)  and  Lilley  (1993).  The 
computed  acoustic  power  agrees  with  the  theoretically  derived  proportionality  to  eMf 
deduced  by  Proudman  (1952);  however  the  constant  of  proportionality  is  smaller  than 
the  analytical  result.  The  acoustic  power  was  computed  from  the  DNS  both  by  ensemble¬ 
averaging  the  instantaneous  acoustic  pressure,  and  also  from  the  two-point,  fourth-order 
correlation  T^^{0)T^x(x).  The  two  methods  for  obtaining  the  acoustic  power  yield  con¬ 
sistent  results. 
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Proudman  (1952)  analytically  obtained  the  following  relation  for  the  acoustic  power, 

Pa  =  a{«V0^ 


where  a  =  13  with  the  particular  choice  /(r)  =  for  the  longitudinal  velocity 

correlation.  In  contrast,  a  ^  2.6  in  the  present  DNS.  The  analysis  of  Lilley  (1993) 
indicates  that  a  is  likely  to  have  a  value  between  3  and  10.  The  analytical  studies  have 
had  to  necessarily  make  some  assumptions  about  the  space-time  statistics  of  the  acoustic 
source  while  the  DNS,  although  free  from  such  assumptions,  is  restricted  to  a  moderate 
to  small  turbulence  Reynolds  number.  The  LES  study  of  Witkowska  and  Juve  (1993) 
which  have  turbulence  Reynolds  numbers  larger  than  those  in  the  DNS  by  an  order  of 
magnitude  give  a  ~  2.5.  Although,  the  DNS  and  LES  studies  give  consistent  results, 
further  DNS  studies  at  higher  Reynolds  number  are  desirable  in  the  future. 
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Fig.  1.  Evolution  of  the 


Fig.  2.  Evolution  of  the 


Fig.  6.  Acoustic  p 
scaled  to  remove 


2 


Fig.  10.  The  frequency  spectrum  of  the  computed  acoustic  pressure. 
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Fig.  11.  DNS  results  at  time 
on  longitudinal  e^ttal  ooxTelationa. 
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Fig.  12.  Effect  of  increased  spatial  and  temporal  resolution  on  F{r\ 
the  longitudinal  spatial  correlation  of  the  acoustic  source. 
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Fig.  13.  Comparison  of  two  methods  of  obtaining  the  scaled  acoustic  power. 
Solid  line  is  an  ensemble  average  from  the  computed  acoustic  pressure, 
and  symbols  are  from  the  computed  two-point  correlation  of  . 
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